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MAGNETIC  SHIELDING  OF  NUCLEI  IN  MOLECULES 
BY  A VARIATIONAL  METHOD* 

by 

James  J?'.  Hornig+  AND  Joseph  O.  Hirschfelder 

Naval  Research  Laboratory,  Department  of  Chemistry 
University  of  Wisconsin,  Iviadison,  Wisconsin 

ABSTRACT 

The  accurate  measurement  of  nuclear  magnetic  moments  requires  a 
calculation  of  the  magnetic  field  at  the  nucleus  arising  from  the  orbital 
motions  of  the  electrons  in  the  external  magnetic  field.  Lamb^  has  calcu- 
lated this  "magnetic  shielding"  effect  for  the  case  of  atoms.  Until  recently 
the  Lamb  theory  has  been  quite  adequate  for  estimating  shielding  effects  in 
molecules,  but  the  accuracy  of  recent  experiments  has  indicated  the  need  for 
a more  comprehensive  theory.  Ramsey^.  3,  4 has  presented  such  a theory 
and  has  discussed  its  need  and  application  quite  thoroughly. 

i 
i 

Ramsey's  results  follow  from  a second  order  perturbation  treatment 
of  the  problem,  and  as  such  are  difficult  to  apply  since  the  wave  functions  of 
excited  electronic  states  of  molecules  are  seldom  available.  In  the  special 
case  of  linear  molecules,  Ramsey  showed  that  the  difficult  terms  could  be 
evaluated  experimentally,  thus  allowing  a determination  of  the  magnetic 
shielding. 

( 

In  this  note  we  outline  a solution  of  the  same  problem  using  the  quantum 
mechanical  variation  principle.  Since  all  results  are  in  terms  of  ground 
state  wave  functions,  such  an  approach  may  facilitate  a purely  theoretical 
evaluation  of  shielding  constants.  As  a variation  function  we  take  a function 

t Z. 

of  the  type  suggested  by  Hylleraas3  and  Hasse  and  applied  successfully  to 
calculation  of  polarizabilities. 


* This  work  was  supported  in  part  by  Contract  N7onr-28511  with  the  Office 
of  Naval  Research. 
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To  calculate  the  magnetic  shielding  at  a particular  nucleus,  we  assume 

with  Ramsey^’  1)  that  we  have  a polyatomic  molecule  in  its  ground 

state,  which  in  the  absence  of  an  external  field  has  no  resultant  electron 

spin  or  orbital  angular  momentum,  2)  all  other  nuclei  have  zero  moment, 

and  3)  the  nucleus  in  question  has  a moment  of  magnitute  JJ.  in  the  same 

direction  as  the  applied  field  H . We  take  this  to  be  the  z direction  and 

hold  the  nuclear  skeleton  fixed  in  an  orientation  specified  by  the  subscript 

A . The  energy  of  the  electronic  system  is  then  calculated  and  all 

/ 

terms  proportional  to  the  product  yuH  are  collected  and  called  WA  . 

The  shielding  constant  is  then  given  by 

w'  v*  H <rx  (i) 

Thus  or^  is  so  defined  that  if  an  external  field  of  unit  strength  is  applied, 

a magnetic  field  with  component  - CT.  parallel  to  the  applied  field  is  in- 

* 7 

duced  at  the  nucleus  by  the  motion  of  the  electrons. 

Using  the  three  assumptions  mentioned  in  the  previous  paragraph, 

8 

the  Hamiltonian  of  the  electronic  system  may  be  written  with  sufficient 
accuracy  in  the  form 


3m 


V*+V  + 5^I 


E(H  + ^)l.i 


(2) 


8mc‘ 


where 
electron, 


ik. 


is  the 


J: 


component  of  the  angular  momentum  of  the  k’th 

J2  f . d . 


7.  For  a more  detailed  discussion  of  the  shielding  constant,  see  Norman 
F.  Ramsey,  Nuclear  Moments  , Wiley  & Sons,  New  York  (1953). 

8.  For  a discussion  of  the  complete  Hamiltonian  of  a system  in  external 
electromagnetic  fields,  see  J.  O.  Hirschfelder , C.  F.  Curtiss,  and 
R.  B.  Bird,  Molecular  Theory  of  Gases  and  Liquids,  Wiley  & Sons, 
New  i'ork  (1954). 
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3. 


Equation  (2)  implies  that  the  vector  potential  depends  only  on  the  external 
field  and  on  the  nuclear  magnetic  moment.  Since  we  will  use  a Hylleraas  - 
Hasse  type  trial  function  and  will  be  interested  only  in  energy  terms  pro- 
px  H > the  parts  of  the  Hamiltonian  proportional  to  jx> 
mjay  be  dropped.  Thus  we  use 


portional  to 
and 


H' 


# ~ + jj.  it'*  H 7t  a H 


where 
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Asa  variation  function  we  take 


where 


"5: 


satisfying 


Then 


is  the  normalized  ground  state  wave  function 

Xt  =E„-£: 
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If  we  choose 


fjrVfr  dc 

<*,/»>  {-$■*  -jr  <j=- 

and  p such  that 

(H)  ° (If)'0 


(3) 
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(5) 


(6) 
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(8) 


and  H 


and  require  that  oC  and  ^ behave  properly  as  jx 
tend  to  zero,  we  will  have  an  approximation  of  the  lowest  energy  level  of 


I 


i i 
\ 


if 

i \ 


i 

! 


II 


Series  3 
WIS-ONR-9 
13  July  1954 

4. 

the  system  in  the  magnetic  field.  For  small  jjl.  and  H the  result 
may  be  expanded  in  a series  in  powers  of  jx  and  H . The  coefficient 
of  the  u.  H term  is  just  the  CT  we  seek. 

v<*,V 

ta.lv  *,]) 


~ 7 2 ' / ^ 

+ Gy  *>» + % •> )( *,  K , #,1  ♦«»[*•»  *■ ] \ 

+ a*,  [<*, 


(9) 


Here  the  symbol 


( ) 


represents  the  matrix  element 


diagonal  in  the  ground  state  evaluated  with  the  nuclear  skeleton  held 
fixed  in  the  orientation  specified  by  “X  , and  ["  0}P^  is  the 

commutator  of  O and  P.  Also, 
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6. 


Thi*  rather  complicated  result  may  be  simplified  somewhat  by  re- 
ferring back  to  our  original  assumption  that  there  was  no  contribution  to 
the  vector  potential  from  the  orbital  motion  of  the  electrons.  Various 

quantum  mechanical  interpretations  of  this  assumption  will  give  varying 

v-  ■ ■ : 

degrees  of  simplification  of  this  result.  The  absence  of  orbital  angular* 
momentum  certainly  implies  that 
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and 


Setting  these  quantities  equal  to  zero  in  affords  a slight  simpli- 

fication. If  vre  assume,  further,  that 

sib  Hz) 
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the  result  becomes  very  much  simpler.  Noting  that  7€t 

exactly,  • 
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(13) 


Consistent  with  Xx  j = 0 , it  may  be  shown  that 
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This  result,  together  with  the  definition  of  the  operator’s  gives,  finally, 
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7. 


Resonance  absorption  and  nuclear  induction  methods  use  molecules 
of  all  orientations,  so  this  shielding  constant  must  be  averaged  over  all 
possible  orientations.  The  first  term,  corresponding  to  the  Lamb  cor- 
rection, is  easily  averaged  since  x,  y,  and  z are  identical  under 
such  an  average.  Finally, 


This  result  may  be  compared  with  the  earlier  results  of  Lamb 
and  Ramsey^.  For  atoms.  Lamb  got 


sL 


Ramsey's  result  for  molecules  is 


a. 


•where  the  sum  over  7^  A extends  over  all  excited  electronic  states. 
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